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1. Introductory. — This paper is a development of two earlier papers, # which for 
brevity I call " Reduction " and " Fitting " respectively. The paper f immediately 
preceding " Fitting" is referred to as " Factorial Moments." 

These earlier papers deal with two problems, which are closely connected and have 
the same solution. For both of them, the data are a set of quantities u , u u u 2 , ... of 
the same kind, which we regard as representing certain true values U , U l9 U 2 , ..., 
with errors e , e x , e 2 , ..., so that u r = U r + e r . These errors may be independent or 
may be correlated in any way. The first problem is based on the assumption (which 
defines the class of cases we are dealing with) that the sequence of U's is fairly 
regular, so that differences after those of a certain order, which we will call/, are 
negligible. This being so, we may alter any u, or any linear compound of the u's, 
such as an interpolation-formula, by adding to it any linear compound of the neg- 
ligible differences. (I use the term " linear compound" in preference to "linear 
function," since there is no consideration of functionality.) The problem is to find 
the value of the resulting sum when, by suitable choice of the coefficients in the 
added portion, the mean square of error of the sum is a minimum. This is the 
problem of " reduction of error." For the second problem it is assumed that U r is a 
polynomial in r of degree j, and the problem is to find the coefficients in this poly- 
nomial by the method of least squares. This is the problem of " fitting." 

The practical solution of these problems for the general case, in which the errors 
are correlated, is not easy. The particular case which is simple is that in which 
the errors all have the same mean square, which by a suitable choice of unit is taken 
to be 1, and the mean products of error are all 0. (In the previous papers I have 
called this system of errors the standard system ; in the present paper the set of u's 
which possesses this property is called a self-conjugate set.) In " Reduction " I gave 
the solution for this particular case in terms of central differences, and in " Fitting " 
I gave the solution in terms of advancing differences and of advancing and central 
sums, formed in a particular way. I also gave expressions in terms of the us, but 
these were rather complicated. It remained to obtain expressions for the mean 
squares of error of the new values, in order to compare them with those of the old 

* " Keduction of Errors by means of Negligible Differences," ' Fifth. International Congress of 
Mathematicians,' Cambridge, 1912, ii., 348-384; "Fitting of Polynomial by Method of Least Squares," 
' Proceedings of the London Mathematical Society/ 2nd series, xiii., 97-108. 

t " Factorial Moments in terms of Sums or Differences," ' Proceedings of the London Mathematical 
Society/ 2nd series, xiii., 81-96, 
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values. In doing this I found that the whole of the work could be very much 
simplified by using certain general theorems, which applied not only to the special 
case of the standard system but also to the general case, and even to a still more 
general problem in which, in the one aspect, the reduction of error is effected by 
means of quantities which are not necessarily a set of differences, or in which, in the 
other aspect, U r is not necessarily a polynomial in r of degree j but is a linear com- 
pound, with coefficients to be determined, of any/+l functions of r ; and the present 
paper is mainly concerned with these general theorems, so that to a certain extent it 
supersedes the previous papers. 

The abbreviations I.e., m.s.e., m.p.e., are used for linear compound, mean square of 
error, mean product of errors. The mean square of error of A is denoted by (A ; A)> 
and the mean product of errors of A and B by {A ;B) or (B ;A). Other special 
notations used in the paper are the same as in the three papers mentioned at the 
beginning of this section, or are explained in §§3, 5 (iii.), 7, 17, and 20. 



Conjugate Sets. 

2. Conjugate Set — (i.) Let A, B, (7, D, ... be a set of quantities, not necessarily all 
of the same kind, containing coexistent errors which are either independent or 
correlated in any way. For the purpose of the following investigations it is 
convenient to consider, in connexion with these quantities, another set of quantities, 
6r, jET, J 9 K, ... , equal to them in number and connected with them by the conditions 
that (l) each quantity of the second set is a I.e. of those of the first set, and (2) the 
m.p.e. of corresponding members of the two sets is 1 and that of members which do 
not correspond is 0. If we replace the quantities of the two sets by A , A u A 2 , ... , 
and 6r , Gf u 6r 2 , ... , we can express this latter condition by saying that m.p.e. of 6? r 
and A s = (s ^ r) or 1 (s = r). The second set of quantities is said to be conjugate 
to the first. 

(ii.) Let the member of the second set which corresponds to C of the first set 
be J. To determine J, let us write 

J = aA + bB + cC+dD+.... 

Then, denoting the m.p.e. of A and B by (A ; B), condition (2) gives 

(A ;A)a + (A ;B)b + (A ; C)c + (A ;D)d+... -0, 

(B ; A)a + (B ;B)b + (B ; C) c + (B ; D) d+ ... -0, 

(C ;A)a + (C ;B)b + (C ; Cf)c + {C ; D) d+... = 1, 

(D;A)a+{D;B)b + (D;'C)c + {D;D)d+...=0 t 

&c. 
2 F 2 
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There are as many equations as there are coefficients a, b, c, d, ... ; and the values 
of these are thus uniquely determined. 

(iii.) The values of a, b, c, ... as found from the above equations have as their 
denominator the determinant 



e 



(B 

(C 



A) {A ;B)(A ; C) (A 
A) (B ;B){B; C) (B 
A) {G ; B) (C ; O) (C 
A) (D ;B){D; C) (B 



D) . . . 
D) . . . 

JL-JT J * • • 

D) . . . 



There is therefore no conjugate set if this determinant is zero. The nature of the 
relations which in this case hold between the errors is considered in Appendix I., § 3. 
(iv.) Since the members of the conjugate set are l.cc. of those of the original set, 
the converse also holds. Regrouping the equations which determine the coefficients, 
it will be seen that the original set is conjugate to the conjugate set ; i.e., that the 
two sets are conjugate to each other. The formulae for the members of the original 
set in terms of those of the conjugate set are 



A = (A ; A)G + (A ; B)H+(A ; C)J+..., 
B = (B;A) G + (B ; B) H+(B ; C)J+... 9 
C = (C;A) G + (0 ; B)H+(C ; C) J+...- 9 

&c. 



"S 



> 



J 



(1) 



These follow from the solution of the equations in (iL), by the ordinary properties 
of determinants; or they may be obtained more simply by determining the 
coefficients of G 9 H, J, ... in each case from the second of the conditions stated 

in (i.). 

(v). By means of these relations we can not only express any I.e. of the quantities of 
either set in terms of those of the other set ? but we can also express any such I.e. in terms 
of particular quantities of one set and those of the conjugate set which correspond 
to the remaining quantities. We can, for instance, express any I.e. of G, H, J, K, ... 
in terms of A, B, J, K, ... by using the first two equations in (l) to determine 
G and H in terms of A 9 B 9 J,K 9 .... The results involve a certain determinant in 
the denominator ; it is shown in Appendix I., § 4, that this is not zero if is not zero. 

(vi.) Two special cases may be mentioned : — 

(a) If the errors of A, B 9 C 9 D 9 ... form a standard system, i.e., if the m.s.e. of 
each of the quantities is 1 and the m.p.e. of each pair of quantities is 0, the conjugate 
set is identical with the original set ; and conversely. A set which is identical with 
the conjugate set will be called a self-conjugate set. 
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(b) If the m.p.e. of each pair of quantities of the original set is 0, but the m.ss.e. 
are not all 1, this is also the case for the conjugate set. The original set being 
A, i?, G, ... ? the quantities of the conjugate set are A /(A ; A), B/(B ; J9), C/(C ; C), ... ; 
and their m.ss.e. are lf(A ; A), l/(B ; 5), l/{C ; C), — 

3. Relations between Original Set and Conjugate Set. — For expressing a member 
of either set in terms of the members of the other set, it is convenient to give them 
a linear order. We therefore denote the members of the original set by S Q9 S l9 S 29 ... S t 
and those of the conjugate set by <r , <r 1? <r 2 , ... o> Also we write 

t r ,t = ni.p.e. of § r and S t = & ir , . . . . . . . . (2) 



n T ,t 



m 



.p.e. of ar r and a- t = fj t$r . . (3) 



(i.) The condition of conjugacy is that (r = 0, 1, 2, ... I; t = 0, 1, 2, ...7) 

m.p.e. of cJ r and <r t = (r ^ t) or 1 (r = t). . . . . . (4) 

(ii.) The expression for S r in terms of the o-'s is (cf. § 2 (iv.)) 



[For, if we write 



then (2) and (4) give 



K = &, 0^0 + fr. 1^1 + fr, 2^2 + • • • + f r, l°"i- • ' • (5) 



5, = a <r + a 1 <r 1 + a 2 (7 2 +...+^o' z , 



f ff t = m.p.e. of ^ and a ar + a l a- l + . . . + a^ 
= «,] 



(iii.) Similarly the expression for cr t in terms of the S'$ is 



o"« = %,A+ , yi,A+%,A+*-*+^A 



(6) 



(iv.) The relations between the f 's and the >/s are easily deduced from (5) and (6). 
If we write 



Z = 



?o,o So.i So, 2 • • • so,2 

M,0 Sl.l M,2 • • • M,J 
S2.0 52,1 S2.2 • • * S2,Z 



M,0 M»l M,2 • • • ?Z,Z 



? V' / 



Z Ptg == cofactor of £ Ptq in Z = Z g tP> . (8) 
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H zz 



»/0,0 %,1 *?0,2 • • • tf<U 



tfl.O ft f l */l,2 • 



^2,0 %,1 % f 2 



• » « 



tfl,l 



*72,J 



W.0 %,1 *7/,2 • * • %.* 



• •••••** 



(9) 



then 



H„ „ = cofactor of 



*lp,q 



m H = H 



g,p? * 



» = Z /Z 

J, — I^mLm ""'"""" ,ri JL * » • • • s • * 



. . . . (10) 

c « * • I JL JL / 

• • • • l A M / 

. . . • [Loj 



(v.) The assumption that there is a conjugate set implies (cf. §2 (hi.)) that 
Z is not = 0. It follows from (13) that H is not = 0. It also follows (see 
Appendix L, § 4 (&)) that none of the principal minors of Z or of H are = 0. 

4. Two Related Pairs of Conjugate Sets.* — (i.) Suppose that there is another set 

of Z + l quantities u 0) u u u 2 , ... u u connected with the 8's by the linear relations 

(r = 0, 1, 2, ... I) 

(14) 



u r = (n) 4+ (n) Si+(r 2 ) 4+ ... + Vi) s t . . 



Then, by the condition of conjugacy of the <f s and the o-'s, 

(r t ) = m.p.e. of u r and ar t . 



(15) 



Let the set conjugate to u 09 u u u 2 , ... u x be y Qi y l9 y 2 , ... y t . Then there are linear 

relations between the y's and the us and between the <rs and the 5s, and therefore also, 

by (14), between the ys and the <rs. To find the cr's in terms of the ys, we write (15) 

in the form 

(r t ) = m.p.e. of o- t and u r ; 

and we see that (t = 0, 1, 2, ... I) 

*t = (o*)y +(i*)yi+(. 2 *)ya+-- + (0y/- ...... (16) 

(ii.) Similarly, if the expression for the 5's in terms of the w's is (s = 0, 1, 2, ... Z) 

S 9 = {s }^o+{si}n i +{s 2 }u 2 + ...+ {s l }u l9 (17) 

where, by (14), 

(r Q )iO t }+(r l ){l t }+... + (r l ){l i } = 0(r 9 *t) or l(r = t) 9 . . . (18) 



then 



W(0,) + {*i}(l f )+...+ {*,}(O = 0(r^0 or l(r = *),. . • (19) 

{sj = m.p.e. of 5 f and y u ....... . (20) 



% = {0«}^o+{l*}o"]+ {2^}(T 3 +... + {Z,}<r,. 



(21) 



REDUCTION OF ERROR BY LINEAR COMPOUNDING. 



205 



(iii.) The above relations can be expressed diagrammatically, thus 



• 


• 


• 


• 


• 


• 


• 
• 


• 
• 


• 


* 


• 




4 


X 


X 


X 


X . . . 


b M 


A 


X 


X 


X 


X • • . 


t? 


*i 


X 


X 


X 


X . . » 


r-l 

b 


*. 


X 


X 


X 


X » . . 


b° 




u 


u l 


u 2 


6vq • ■ • 






o 


T-i 


<M 


CO 






^ 


55i 


Sfc 


^i ... 





The crosses represent the ( ) coefficients if they are the coefficients of the S's in 
the us and of the y& in the o-'s, and the { } coefficients in the converse case, 
(iv.) Similarly, if we write (r = 0, 1, 2, ... I ; t = 0, 1, 2, ... I) 



then 



{ft] = na.p.e. of y r and ar t , 



Vr = Oo] ^o + M 4 + W 4 + • • • + |>J 4 
<"■« = [°J M o+ [1 J «i + [2,] w 2 + . . . + [ZJ u, 



(22) 

(23) 
(24) 



5. Sums as Conjugates of Differences. — The cases of importance are those in 
which the <?s are successive differences of the u'&. It will be found that in these 
cases the o-'s are l.cc. of successive sums of the g/s. 

(i.) Let the 5's be the advancing differences of the us, i.e. 

Oq == wqj 0| == ZA Uq) • . . O r == l\ Ufa • • • . 

Then the diagram for the ( ) coefficients is 



• 
• 


• 
• 


• 
• 


• 
• 


• 
• 


• 
• 


• 

3» 


• 




• 




• 




• 

1 

-JL « • • 


• 

CO 

b 


4 








1 


O • • * 


£ 


A . 





1 


2 


O ". • • 


b 


So 


i 


1 


1 


-*- • ft,, • 


b° 




u a 


u x 


u 2 


wg • • • 






o 


r-l 


w 


CO 






5Si 


S3> 


5* 


^i • • • 
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so that 



and, generally, 



or 



f 



*0 = ^0 + 2/1 + 2/2 + 2/3+-. -+Vl> 

a-! = y 1 + 2y 2 +Sy B +4y A +...+ly l9 

0-2 = y2 + 3y z +6y l +l0y 6 +...+$l(l-l)y l9 . 

= y / +(/+l,l)y /+1 +(/+2,2).^ +3 +... + (Z,Zr-/)^ 
= 2 (q,f)y q ; ............ 



or, in the notation of "(Fitting," § 4, and "Factorial Moments," 



oy = S^V/. 



• * 



(25) 



• • • ft • ft • I j£J \J / 



The cms can be obtained by successive summations of the ys in reverse order, i.e. 
from y t to y d , as shown in the following diagram, in which the crosses represent 
entries in a sum- or difference-table : — • 











Vl 





yi 

x 

x y^ 

Ox x 







Vl-2 



x 



X 



<*h 



X 



2/ 2 



or a 



0s 



2/i 



OTo X 



2 



0*1 



2/0 



* #0 



or, 







U 



■0 = ^0 



*1 



1 u 2 

x 4 

^o X (I 



X 



u 



1-2 



u 



z-1 



Ut 



X 



*l-l 



*, 



(ii.) Let the S's be the central differences of the us. Then it will be found in the 
same way that 

(a) If the us are u 09 u l9 u 2 , ... u 2n , so that 
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then 



<*2h 



r — n 

2L 

r= —n 



r 



%h)y n 



+r> 



(9*7\ 
* • • • • • • • \ * / 



r = n 



<J 2h-l 



2 (r,2h-l]y n+r ; (28) 



r — —n 



(b) If the u 9 & are u m u l3 u 29 ... u 2n __ l9 so that 



then 



S = y.u n _& S t — Su n _ b $ 2 = v.8 2 u n _ h S 3 == Shi n _ b ... 9 



r — n 



°2A-1 



2 [r- 



o , ^jA 



r = -?i-f 1 



r = n 



•*•/ t/n+r—1* 



<t-m-2= 2 (r-i,2A-2]y B+r _ 1 . 



fSK — n+1 



. . . (29) 



(30) 



(iii.) The values given by (25)-(30) may be expressed in terms of successive sums 
by the formulae given in "Factorial Moments." The notation, however, can be 
simplified. Suppose that we have a set of quantities . . . F 09 F l9 F 29 ... corresponding 
to values ... 0, 1, 2, ... of some variable, and that we form the table of successive 
differences (and also, if we like, of successive sums) of the F's. Then the Lagrangian 
formula for F in terms of F p9 F p+l9 ... F p+t9 which can be expressed in a good many 
different ways, may be regarded as the formula for it in terms of the whole 
(unlimited) set of differences (and sums) which form a triangle with its apex 
at A*F P ; and we can denote it by L {F e ; 6?F P }. With this notation, the above 
results may be written 



(25) 
(27) 
(28) 
(29) 

(30) 



o> — 



7 



x{(-)/2^V/;2 % }];:^ (si) 



1 -itszn+h 



cr 2h = [L {v°- 2h+l y n ; <ry n+t }j t= :; n ^ 



°"2A-i — L-^ \—<r y n ; <ry n+t }] t;s: __ n __p 



. . 



. . . i iojuj 



(33) 



<rih-i = [L {—fx(r ih y n ^ ; ary n+t _ 1 }'] 



— fn _L A 



t = n+ 



"71 + 2 



A? 



°"2h-~2 — 



L {<r 2h Vn-| I Vl/n + t-l}. 



t-n + % 



. . . . (34) 
. . . . ^oo j 



The L is distributive as regards the first member inside the { } ; e.g., in the case 
of (31), 

A<r a +B<r, = [L {A2»y a -B2% ; 2y ( }]' 



t = ' 



(iv.) More generally, suppose that the S'b are the successive differences of the us 
according to any system of differences ; by which we mean that S 8 is either a definite 
difference of the us of order s (the us themselves being differences of order 0) or a 
I.e. of such differences. Then (r t ) of ( 14) is a polynomial in r of degree t 9 and a- t is 

vol. ccxxi. — a. 2 G 
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of the form 2 <f> t (q) y q9 where <j> t (q) is some polynomial in q of degree t. It follows 

g = 

that any I.e. of o- , <r u cr 29 ... cr t is also of this form. 

(v.) If we denote the f th moment of the ys by Mp then M f is of the form 

2 fy (q) y q . Hence a- t is a I.e. of M 09 M l9 M 29 ... M t ; and M t is a I.e. of a- Q9 o-j, <x 2 , ... a> 

More generally, any I.e. of <r , a- 1? a- 2? ... cr^ is a I.e. of M 09 M l9 M 29 ... M t ; and 
conversely. 

Kepuction of Error (General). 

6. General Theorems. — Let A 9 P, G 9 ... P, Q 9 R 9 ... be a set of quantities as in § 2, 
but all of the same kind. If 

w = a^4 + &P + c(7+... (with or without terms in P, $, P, ...), 
cc = w+jpP+gQ+ri2+... , 

where a 9 b 9 c 9 ... are fixed and p 9 q 9 r, ... are arbitrary, and if we choose p 9 q, r, ... 
so as to make the m.s.e. of as a minimum, the resulting value of x is called the 
improved value of w, using P, Q, R, ... as auxiliaries. The following are general 
theorems; some are quite elementary, but it is convenient to state them here. The 
specially important theorems are (III.) and (XIII. ). The assumption mentioned 
under (VI.) should be noted. If strict proofs of (I.) and (II.) are required, the 
method should be that of Appendix I., § 2. 

(I.) The m.p.e. of A and any I.e. of A, B, 0, ... is the same I.e. of the m.pp.e. of 
A and A, B, C, ... [i.e., m.p.e, of A and aA + bB + cC+ ... is a (A ; A)+b(A ; P) + 

c(A ; G) + ...]. 

(II.) The m.s.e. of any I.e. of A, B, 0, ..., or the m.p.e. of any two such l.cc. y is 

found by squaring the former or multiplying the latter and replacing squares and 

products by the corresponding m.ss.e. and m.pp.e. [i.e. 9 m.s.e. of aA-\-bB + cC+ ... 

= a 2 (A ;A)'+2ab(A ; B) + & 2 (B ; B) + 2ac(A ;G)+2bc(B; G)+c 2 (C; (7) + ..., and 

similarly for m.p.e. of aA + bB + eC+... and a'A + b'B + c'G+ ...~]. 

(III.) If the improved value of A, using certain auxiliaries, is A + a, then the 
m.p.e. of A + a and eacA of the auxiliaries or a or any other I.e. of the auxiliaries is 
zero. [Let the auxiliaries be P, Q 9 P, ..., and let J. + a = A+pP + qQ + rR+ ... . 
Then the m.s.e. of A+(p + 6) P + qQ + rR+ ...(= A + a + 6P) is (^4 + a;J+a) + 
29 (A -{-a ; P) + 2 (P ; P). In order that this may be a minimum for = 0, 
(A + a ; P) must be zero. Similarly for (A + a ; Q) 9 (A + a ; P), ... . This proves the 
first part of the theorem ; the second then follows from (I.).] Hence 

(IV.) i/ fcAe improved values of A erne? o/ B, imwgr m eacA case £Ae scwe set of 
auxiliaries, are A + a and B + /3, ^e m.pp.e. of A + a and B + /3, o/ A + a and B, and 
o/ A and B + /3, are aZ£ equal; and 
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(V.) If the improved value of K, using certain auxiliaries , is A + a, and that o/*B, 
using some only of these, is B + 8, the m.p.e. of A + a and B + /3 is equal to that of 
A + a and B. 

(VI.) If the improved value of A, using P, Q, R, ... as auxiliaries, is 
A + pP + qQ+rR+... , the values of p, q, r, ... are given by a set of linear relations 
between p, q, r, ... and the m>.pp.e. of A and P, Q, R, ... . [The equations are given 
by (III.), viz., since (A+pP-t-qQ+rR + ... ; P) = 0, &c., 



(A ; P) +p (P ; P) +q(Q ; P) + r (P ; P) + ... 
(A ; Q)+p(P ; Q)+q{Q ; Q)+r(R ; Q) + ... 
(A ; R) +p (P ; R) +q(Q ; R) + r (P ; P) + ... 

&c. 



0, 

o, 



It is assumed that the determinant 



(P ;P) (Q; P) (R;P) . . . 

(P ; Q) (Q ; Q) (R ; Q) . . . 

(P ; i2) (<? ;R) (R;R) . . . 



is not zero ; i.e. (see § 2 (hi.)) that there is a set conjugate to P, Q, R, ... .] 

(VII.) jpbr any particular set of auxiliaries there is one and only one improved 
value of A. [This follows from the fact that the equations in (VI.), on the 
assumption there stated, have one and only one solution.] Hence we get the 
converse of (III.) : — - 

(VIII.) If x is the sum of w and a Lc. of P, Q, R, ... , and if the m.p.e. of x and 
each of P, Q, R, . . . is zero, then x is the improved value of w, using P, Q, R, . . . as 
auxiliaries. As a particular case : — 

(IX.) If the m.p.e. of w and each of P, Q, R, ... is zero, the improved value of w, 
using P, Q, R, ... as auxiliaries, is the same as the original value. 

(X.) The improved value of P, using P, Q, R, ... as auxiliaries, is P— P = 0. 
[This follows either from (VIII.) or from taking A = P in the equations in (VI.).] 
Hence 

(XL) If w is a I.e. of A, B, C, ... P, Q, R, ... , the improved value of w, using 
P, Q, R, ... as auxiliaries, is the same as that of the quantity obtained by adding to 
w any I.e. of P, Q, R, ... ; and, conversely : — • 

(XII.) If the improved values of w and of w', using in each case the same set of 
auxiliaries, are identical, then w and w' either are identical or differ by a I.e. of the 
auxiliaries. 

(XIII.) If the improved values of A, B, C, ... , using in each case the same set of 

2 G 2 
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auxiliaries, are A -ha, B + /3, C + y, ... , the improved value of any I.e. of A, B, C, ... , 
using these auxiliaries, is the same I.e. of A + a, B + /3, C + y, ... . [Let the I.e. of 
A, B, C, ... be w = a^4 + 6jB + cC y + ... . We want to prove that x = a(A + a) + 
b (P + /3) + c(C-hy) +... is its improved value. We can do this in either of two 
ways : 

(i.) By (III.), the m.p.e. of x and each of the auxiliaries P, Q,„ R, ... is zero ; and x 
differs from w by a I.e. of P, (?, R, .... Hence, by (VIII.), as is the improved value 
of w, using P, Q, R, ... as auxiliaries. 

(ii.) A more direct proof follows from the linearity of the equations mentioned in 
(VI.). It is not necessary to set out the proof here.] 

(XIV.) If A, B, ... 0, D, ... P, Q, ... R, S, ...fall into two classes A, B, ... P, Q, ... 
and C, D, ... , R, S, ... , such that the m.p.e. of each member of the one class and 
each member of the other class is zero, then the improved value of a I.e. of any of the 
members, using P, Q, ... R, S, ... as auxiliaries, is to be found by taking the two 
classes separately, i.e., by using P, Q, ... as auxiliaries for the terms in A, B, . . . P, Q, . . ., 
and R, S, ... as auxiliaries for the terms in C, D, ... R, S, ... . [For the m.s.e. of 
aA + bB+ ...'+cC+dD+ ... + pP + qQ + ... + rR + sS+... is the sum of those of 
aA + bB+...+pP + qQ + ... and cC+dD+ ...+rP + s>S+ ... , since the m.p.e. of these 
latter is zero; we cannot reduce the m.s.e. of the first of them by adding terms in 
R, S, ... , or that of the second of them by adding terms in P, Q, . . . ; and the result 
is therefore the same as if we considered them separately.] 

(XV.) If the improved value of w, using P, Q, R, ... as auxiliaries, is 
x = w + pP + terms in Q, R, ..., this is also the improved value of w + pP, using 
Q, R, ... as auxiliaries. [For x differs from w+pP by terms in Q, R, ..., and the 
m.p.e. of x and each of Q, R, ... is zero.] This can be stated more generally as 
follows : — 

(XVI.) If the improved value of A, using a set of auxiliaries S, is A + a, and if 
we divide S into two sets, S x and S 2 , and the corresponding parts of a are <Xj and a 2 , 
then A -fa is the improved value of A + a 1? using S 2 as auxiliaries. [We may take 
... P, Q to be S l9 and R, ... to be S 2 . The theorem states that, if the improved value 
of A, using ... P, Q, R, ..., is A+ ... +pP + qQ + rR+ ..., this is also the improved 
value of A+ ... -i-pP + qQ, using R, ....] 

(XVII.) The following corollaries of (III.) may be noticed, though we shall not 
require them. If the improved values of A and of B, using in each case the same 
set of auxiliaries, are A + a and B + /3, then 

(1) (A-ha; A+ol) = (A ; A)-(a; a) 
and 

(2) (A + a ; P + /3) = (A; B)-(a; 0). 

.7. Notation: and Particular Values. — (i.) It will now be convenient to adopt a 
linear arrangement of the quantities we are dealing with, and we therefore replace 
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A, B 9 G 9 ... P 9 Q 9 R 9 ... by S 09 S l9 S 2 , ... § j+l9 § j+29 ... § t . The order is quite arbitrary, so 
far as any general theorems are concerned ; but it will usually be convenient to place 
the auxiliaries last. If, for instance, we are using all but j+1 as auxiliaries, we 
denote those not so used by S 09 S l9 $ 29 ... S j9 and the auxiliaries by ^ +1 , S j+2 , ... S t ; the 
improved values are then denoted by ( ).. 
(ii.) We use the following notation ' — 

(e^j = improved value of Sy, using S's after Sj ; 
E ; = (e^ = improved value of S j9 using all subsequent <^s ; 
(X /t ^ = m.p.e. of (e f )j and (e g )j ; 

Aj = (Xj, j)j = m.s.e. of E^. 

Where there is no doubt as to the <5s used as auxiliaries, (e f )j and (\y t9 )j can be 

replaced by e f and Xy >r 

(iii.) By(X.)of§6— 

(e / ) j = 0if/>i; (36) 

(X /f ^=0.if/>yorflr>y (37) 

(iv.) By(IV.)- 

( X /.*)j = m 'P' e - of ( € f)j and ^ 1 . / 3g x 

= m.p.e. of Sy and (e g )j J 
A^ = m.p.e. of E, and Sj . (39) 

8. Improved Values in terms of Conjugates.- — In "Fitting" I have given some 
formulae for improved values in terms of sums. These may be regarded as derivable 
from a general theorem relating to the expression of improved values in terms of 
members of the conjugate set. The theorem is given by (XIX.) and (XX.) below; 
(XVIII.) is a particular case. 

(i.) Take any one of the <$'s as S . By (6) — 

Hence <r fa 0t0 differs from S by a Lc. of the other ci's. Also the m.p.e. of <t / %)0 and 
each of these other 8's is zero. It follows from (VIII. ) of § 6 that o- fa 0t is the 
improved value of S 09 using the other <?s as auxiliaries. The m.s.e. of this improved 
value is *7 ,o/U,o) 2 = l/% - Hence— 

(XVIII.) The improved value of any member of the original set, taking all the 
other members as auxiliaries, is the product of the corresponding member of the 
conjugate set by a constant ; this constant being — the m.s.e. of the improved value. 

(ii.) The first part of the more general theorem is :■ — • 

(XIX.) The improved value of any I.e. of a set of quantities, using those after the 
first j + 1 as auxiliaries, is a I.e. of the first j+ 1 of the conjugate set. 
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For, if w is a I.e. of S , S^ § 2 , ... S l9 we can (see §2 (v.)) express it as a I.e. of 
<*b>"°"i> 02, ••• <Tj 9 ■ Sj+i, 8j+to'-'$i' Let the result be (2) + (A), where (2) is a I.e. 
of <r , cr u o- 2 , ... o^, and (A) is a I.e. of S j+1 , § j+2 , ... § t Then (2) differs from w by a 
I.e. of these latter S's, and the m.p.e. of (2) and each of these 8s is ; hence, by 
(VIII ), (2) is the improved value of w, using these 8s as auxiliaries. 

(hi.) Further — 

(XX.) The coefficients of the o-'s in the improved value of the I.e. are the m.pp.e. of 
this improved value and the corresponding 8 s. 

For, if the improved value of w is x, and we write 

X = b cr + h l (T l + b 2 cr 2 4- ... ^rbjCTj, 

then, by the condition of conjugacy of the <r s and the 8s, 

m.p.e. of x and § f = &/. 

(iv.) This would give us a solution of the problem of finding the improved value, 
if we could find the m.pp.e. Ordinarily, w is or can be expressed in terms of the cf's, 
and we do not find its improved value independently, but deduce it from those of the 
8s up to Sj. The improved value of 4 is, by (iii), 

( € h)j = (\h)j°-o + {Kh)jVi+ ••• +(\,h)j< r j> (40) 

and the m.pp.e. that we really require are therefore the values of (X/, A )j. With 
regard to this, see §9. 

(v.) As the converse of (XIX.) it may be noted that — 

(XXI.) A quantity of the conjugate set, or a I.e. of such quaritities, cannot be 
improved by means of the non-corresponding quantities of the original set; e.g., 
a I.e. of cr 3 and <r 4 cannot be improved by using the 8s, other than S 3 and S 49 as 
auxiliaries. This follows from (IX.) of §6, since the m.p.e. of S r and <r s is 
unless r = s,. 

(vi.) If, as in § 5, there are related conjugate sets of us and ys, and the 8s are the 
differences of the us, it follows from § 5 (v.) that (2) in (ii.) above is a I.e. of the 
moments of the ys up to the j th . (XIX.) is therefore a generalisation of the 
theorem, for a self -conjugate set, that the improved values are l.cc. of the moments ; 
and, in fact, it explains the appearance of the moments in this connexion. 

9. Mean Products of Error of Improved Values. — (i.) We have found, in § 8 (iv.), 
that 

( € h)j = (\ h)j °"0 + (\ h)j <Tl + ' ' • + (\", h)j <*> 

To obtain the A's, we introduce the condition that this shall differ from S h by a I.e. 
of 8 b after §j. 
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(ii.) Substituting for cr , <r l9 <r 2 , ... from (6), this condition gives 

+ (*i, h)j (%, A + m, A + • • • + nu A) 
+ ... 

= 4+ terms in 5 i+1 , 5 i+3 , ... 4 
Equating the coefficients of '5 , <^, 5 2 , ... ^, we. find that (/= 0, 1, 2, ...y) 



V. 



;o(Ao,^ + ^i(Ai,A-+... + ^i(^^- = 0(/^^)or l(/= A). . . . (41) 



Let us write 



H, = 



>?0,0 ^0,1 ^0,2 • ■• • *l0 t j 



m t o ^1,1 ?i,a • 

%,0 %,1 ^2,2 • 



fl.j 



%.j 



fj.o ^m fya 



7;. 



y.j 



H 



£> h t j 



Then 



*/. 



and therefore, by (41). 



\^,a)j — H^, ^ ; /H ; . 



• • . • 



(42) 



cofactor of ^ in Hj = H MJ . (43) 



;oH 0)M + ^ 1 H 1)M +...-f^H j} ;, j = 0(/V A) or K 3 (f= h) ; . . . (44) 



(45) 



Substituting in (40), we obtain (e k )j in terms of the o-'s. 

For the particular case of g = h =j 9 (\, k )j becomes A p and H^ tM becomes H^ ; 

so that 

A j = HjjRj • (46) 

(iii.) As an example, suppose that we have several independent observations, of 
unequal accuracy, of a single quantity U, and that we wish to obtain a suitably 
weighted mean, which may be regarded as the improved value of any one of the 
observations. Let the observed values be 'u > u u u 2 > ... , and their m.ss.e. 
a 2 9 af, a 2 2 , ... ; the m.pp.e. being 0, since the observations are independent. We 
take § to be one of the us, and $ u § 2 , ... to be its successive differences. Then j = 0, 
since the true values of the first and later differences are all 0. Hence, by § 8 (L), 
the improved value is cr /(m.s.e. of o- ). But, by § 5 (i.), cr = 2y ; and, by § 2 (vi.) (6), 
y r = u r /a r \ so that m.s.e. of <r = 2l/a 2 , Hence the improved value is 
2(w/a a )/2(l/a a ). 



214 DE. W. F. SHEPPAED ON 

(iv.) When j is relatively large, the solution given in (ii.) above can only be 
regarded as a formal one, since it involves calculation of determinants. I have not 
been able to provide a general solution which shall avoid determinants ; but it will 
be seen in §§ 17-19 that, if we can find the values of certain quantities occurring in 
the formula, we can deduce the X ? s and thence the coefficients of the cr's. These 
latter are important as giving us formulae which contain only a few terms and are 
therefore suited for numerical calculation. 

10. Expressions in terms of a Related Set.— Suppose that there is another set of 
Z+l quantities u 09 u l9 u 29 ... u h connected with the 5's by linear relations ; and let the 
set conjugate to the us be y 09 y l9 y 29 ... y z . We shall take the relations between the 
us and the (?s and between the ys and the 8's to be, as in § 4, (r = 0, 1, 2, ... I) 

u r = (n) s + (n) s x + (r 2 ) s 2 +... + (n) s l9 ...... - (47) 

», = W»o+Wi+Wi J +... + Wi,. , ..... (48) 

(i.) Let the improved value of u r , using $'s after S j9 be v r . Then, from (47), by 
(XIII. ), remembering that, by (36), 

(^ = oif/>y, 

we have 

^==(n)(4- + ( r i)(4- + --- + ( r 0(4 • • ( 49 ) 

= (n) (<*); + (n) (ejj + ...+ (rj) (ej)j (49A) 

Thus the vs are related to the e's in the same way that the us are related to 
the S's. 

(ii.) Similarly, if the improved value of y r9 using <5's after S j9 is z r9 we have 



z. 



= Oo] (*o)j + Cn] (ei)j + . . . + M (4 (50) 

= W («o); + W (4- + • • • + W (4 ; (50a) 

and the zb are related to the e's in the same way that the ys are related to the S's. 

(iii.) Let w be any l.c. of the 8s or of the us or ys 9 and let x be its improved 
value, using S's after Sj. Suppose that x is expressed in terms of the ys 9 the 
coefficients being p 09 p l9 p 29 . . . p l9 so that 

v=p*y*+Piyi+p& 2 +...+p l y l \ ...... . (51) 

and let 

X^ = m.p.e. of x and (e g ) j9 
so that 

X^ = if g > j. 

Then, by the condition of conjugacy of the us and the ys, 

m.p.e. of x and u r = p r . (52) 
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Hence, by (IV.) and (49) or (49a), 

p r = m.p.e. of x and v r 

= (r )*o + (r l )\ l + (r 2 )\ 2 +...+(r l )\ l (53) 

= (n)\) + (n)A 1 + (r 2 )\ 2 + ... + {r j )\ j . . (53a) 

Thus the ps are related to the A's in the same way that the vs are related to the es, 
or the u's to the <5's. 

(iv.) Similarly, if 

x = q u + q 1 u l + q 2 u 2 +...+ q x u h 

the qs are related to the A/s in the same way that the zs are related to the e s, or the 
ys to the <5's. 

11. Special Case of Differences. — The important practical case is that in which the 
$& are successive differences of the us> in the general sense explained in §5 (iv.). If 
the differences of order exceeding j are negligible, we can use them as auxiliaries for 
improving the us or the <5's or any l.c. of the us or the cf's. 

(i.) Since the $& are successive differences of the us, (r t ) is (§5 (iv.)) a polynomial 
in r of degree t. 

(ii.) By § 10 (i.) the e's are the differences of the vs according to the same system ; 
and v r is a polynomial of degree j in r, the differences of the vs of order exceeding j 
being zero. 

(hi.) With the notation of §10 (hi.), the X's are the differences of the p's according 
to the same system ; and p T is a polynomial of degree j in r, the differences of the 
ps of order 'exceeding j being zero. 

(iv.) If we form the differences of the us in the usual way, there will be I differences 
of order 1 , 1 — 1 of order 2, and so on. The l—j+1 of order j\ namely A j u , A j u 2 , . . . A j Ui_ j9 
will differ from one another by l.cc. of the differences of higher order ; and therefore, 
by (XL), they will have the same improved value. If w~e denote this by E, then, 
if w = pA%o + gA% 1 + rA% 2 + ..., the improved value of w is (p + q + r + ...)E. 

Relation of "Reduction of Error" to "Fitting" (of a Polynomial). 

12. Standard System. — In the case of a standard system, the process of reduction 

of error and the process of fitting a polynomial (by least squares or by moments) give 

the same result. The following is a proof of this, not involving the properties of 

conjugate sets. The observed values are taken to be u , Ui, u 2y ...u t ; and 2 denotes 

t 
summation for t = 0, 1, 2, ... I. 

(i.) If the polynomial which we are fitting to the us is 

v q — a + aiq + a 2 q 2 + ... H-o^', ......... (54) 

VOL. CCXXI. — A, 2 H 
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the values of the as when we fit by least squares are given (" Fitting," §§1, 2) by the 
equations (f= 0, 1, 2, ... j) 

2g^. a + Xqf + \ a x + ... + 2g /+y . a j = 2g%, = M, (55) 

These are the same equations that are given by the method of moments. 

(ii.) The above equation (55) is a statement that the f th moment of the vs is equal 
to that of the us. In order to prove that the process of reduction of error, using 
differences of order exceeding/ as auxiliaries, gives the same result, it is sufficient to 
show (a) that the improved value of u q as given by this process is of the form of v q 
in (54), and (0) that the f th moment of the improved values of the us is equal to 
that of the original values for/ = 0, 1, 2, ...j. 

(hi.) We have shown, in §11 (ii.), that the improved value of u q is a polynomial of 
degree y in q. This establishes (a). 

(iv.) By (XIII), the/ th moment of the improved values of the us is equal to the 
improved value of their f th moment. In order to show that this is equal to the 
original value of the f th moment, it is sufficient, by (IX.), to show that the m.p.e. of 
the original / fch moment and every difference of order exceeding j is zero. 

(v.) Let the k th difference of u r _ k be 

S k = Jc u r — k#i rmmX + ... 4- ( — ) k \u r „ k . 
Then the^ 11 moment is 

... + r f u r + (r— l) f u rMl + (r— 2) f u r _ 2 + ... , 
and the m.p.e. of this and S k is 

k r f -k i (r~l) f i-k 2 (r~-2y- ... . 

But this is the k tK difference of (r—k) f , and is = if k >/. 

This proves the proposition. 

13. Fitting by Least Squares.— Next suppose that the set is not self-conjugate. 
If the cfs were the differences of a set of us, we should fit a polynomial of degree 
(say) y to the u's. This suggests that, in the more general case, the us being 
connected with the 8's, as in § 10, by the relation (r = 0, 1, 2, ... I) 

w r = (r )a + (r 1 )i 1 + (r a )*, + ...+(r l )^ (56) 

we should try to fit an expression of the form. 

v r = (r )e ( > + (r 1 )e 1 + (r 2 )e 2 +..'.+(r J )ej (57) 

to the us by an appropriate method of least squares ; the (r)'s being the same as in 
(56), and the e's being the quantities to be determined. 
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(i.) If the ys are conjugate to the u% and if 

Y/y, s ~m.p.e. of y T and y s = ^ r , . (58) 

then (see Appendix II.) the direct (or a priori) probability of the occurrence of the 
given set of it's, if the vs as given by (57) were the true Z7 7 s, is proportional to 

exp {-j22^ M (w f -v r ) K-^)}, 



r s 



where 2 denotes summation for t = 0, 1, 2, ... I. The principle of the method of least 
squares therefore leads us # to choose the e's so as to make 



r s 



a minimum. Differentiating with regard to each of the es, this gives (f = ? 1,2, . . . /) 

2 {(°/)^o.. + (l/)V r i.«+--- + (^)V r /..} ( v s ~~ u s) = (59) 

s 

But, by (58) and (16), 

(°/) ^o,, + ( x /) V r i,.+ ••• + ('/) ^.. = m -P- e - of V, and (0/) 2/o + (l/) &+.-. + (£/) ^ 

= m.p.e. of ^ s and oy (60) 

Denoting this, as in § 4 (iv.), by [fy], the equations (59) become '(/ = 0> 1> 2, ...jT) 

2[s/](^—^) == 0. (61) 



(ii.) Instead of fitting an expression of the form given by (57) to the us we might 
fit a corresponding expression to the ys. Since 



v* = W ^o+ W ^i+M 4+ ... + W a 

the expression to be fitted would be of the form 



* a « » 



■ (62) 



«, = [*o] e + [*J *i + [* J e 2 + . . . + [* 



6j. . 



• « • . . \\Jtjt 



* Strictly, we ought to choose the e's so as to make B . (7 exp - |P a maximum : where 

P is the direct probability of occurrence of the particular e's denoted by € , «i, € 2 , . . . ej, and is therefore 
some function of these latter ; and C exp - \P is the direct probability of occurrence of the particular 
values of u -~ v on the assumption that these values of the e's are the correct ones, C being some function 
of these c's. But I have assumed, as is commonly done, that the range of practically possible values of the 
e's is so small that B and G may be treated as constants, so that we have only to consider the maximum 
value of exp ~ |P. 

2 H 2 
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We should have to choose the e's so as to make 

2Sff r , t (y r -a r )(y,-2,) 

r s 

a minimum, where 

?r n? EE m.p.e. of u r and u s (64) 

This would give 

2(r / )(« r -y r ) = (65) 



r 



(iii.) The e's given by (65) are the same as are given by (61). For we have seen in 

§ 4 that 

(r f ) = m.p.e. of u r and ay, (66) 

[sy] = m.p.e. of ^ and ay (67) 

If we express the uh in terms of the $'s, and write 

8 t 

then, by the condition of conjugacy of the o-'s and the 8's, and by (66), 

A t = m.p.e. of cr t and 2 [ay] u s 

s 

= 2 [*,] (a,). 

This is symmetrical, and we should get the same expression for the coefficient of S t 
in 2 (r f ) y ri so that 

2 [^u. = 2 (ry) y r (68) 



r 



Similarly, if we substitute the values of v s from (57) and of z r from (63) in 2[sy] v s 

s 

and in 2 (r f ) z r9 the coefficients of e t in the resulting expressions are equal. Hence 

r 

(61) and (65) are identical equations in the e's. 

(iv.) The values of the e's as given by these equations are in fact independent of 
the u'& or the y'a. For the value of A t as found in (iii.) above is 

2 [_Sf] (s t ) = m.p.e. of 2 [fy] u s and 2 (s t ) y s 

S S 8 

= m.p.e. of ay and <r t9 ' (69) 

by (24) and (16). Hence, denoting the m.p.e. of ay and <r t9 as in § 3, by rj ftt , the e's 
given by (61) or (65) are the same as would be given by (/ = 0, 1, 2, .../) 

2 *lf t t € t ~ 2 rjf t o t (7 0) 

t-o t=o 
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(v.) The ordinary method of least squares would consist of making 2 (v 8 —u s ) 2 a 

s 

minimum, and would lead to equations 

2(s/)(^-<) = 0, 

s 

which would not give the most probable values of the e's. ' 

14. Fitting by Moments. — (i.) The ordinary method of moments, adapted to the 
case in which the S's are not necessarily the successive differences of the us, would 
consist in equating the values of 2 (s f ) v s and 2 (s f ) %i s . This, as will be seen from 

§13 (v.), would not give the most probable values of the es. 

(ii.) In order to obtain the most probable values of the es by equating moments of 
the vs and of the it's, we must write (say) 

M f = X[s f ]u s , . . (71) 

and define the f th moment of the us as being M f or a definite I.e. of M f , M f _ u 
M f __ 2 , ...M . But the coefficient of u s in M f would then not be given definitely by 
the relations between the us and the <^'s, but would depend also on the law of 
correlation of errors of the us. We see, however, from § 13 (hi.), that we have also 

M f = S{r f )y r (72) 



r 



and that we get the same result by equating moments, defined in this way, of the ys 
and the z'&. In the ordinary case in which the <5's are successive differences of the 
us, the coefficients of the ys in (72) are binomial coefficients, and the ordinary 
moments fall within the definition given above. It follows that in fitting a 
polynomial to a set of quantities (not being a self-conjugate set) by the method of 
moments, the moments which ought to be equated are not those of the quantities 
themselves and their assumed values, but those of the conjugate set of the former and 
the corresponding l.cc. of the latter. 

15. Reduction of Error. — If we improve the S's or the us by means of the §'s 
after S j9 the improved values of these latter are zero, and those of the § 7 s up to Sj are 
obtainable from^XXI.) of § 8, which states that the improved values of the o-'s from 
or Q to o-j are equal to the original values. Using (6), this gives (/= 0, 1, 2, ... j) 

t-j t = i 

2 %t( € t)j= ^ f] ftt S t (73) 

Comparing this with (70), we see that the es given by this process are the same as 
those given by the process of fitting the expression in (57). 

16. Difference of the Tivo Processes. — Although the two processes lead to the 
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same result, they are essentially different. This is explained in § 22 of " Reduction." 
The main difference may be expressed as follows : — 

(i.) In "fitting" we deal directly with the particular case. We assume that the 
true values follow a specified law, involving unknown constants, and we deduce 
values for these constants from the data by the principle of inverse probability. 

(ii.) In "reduction of error" we do not use inverse probability, and we only deal 
incidentally with the particular case. We regard the aggregate of the data as one of 
an indefinitely great number of possible aggregates from the same true values, and 
we use a method which, will reduce as much as possible the m.s.e. of these possible 
aggregates. 

Some Steps in the General Solution. 

17. Preliminary. — (i.) Our object is to find the improved value of any I.e. of the 
(fs or the us, and the m.s.e. of this improved value or the m.p.e. of two improved 
values. Ordinarily, as already stated in § 8 (iv.), the quantity to be improved would 
be expressible in terms of the cf's, so that we need consider only the improved values 
of the cf's, i.e., the Vs. There are then four problems before us, viz. : (l) expression 
of the e's in terms of the <^'s ; (2) expression in terms of the o-'s ; (3) expression in 
terms of the ys ; (4) determination of the A's. For practical purposes (2) is more 
important than (l) or (3), since there will be fewer coefficients involved. 

(ii.) Although it does not seem possible to obtain a general solution, otherwise 
than by determinants, there are some general propositions that indicate stages in the 
solution. If, without necessarily finding the complete expressions of the e's in terms 
of the cJ's, we can find for each e the coefficient of the first of the auxiliaries, then it 
will be seen from § 18 that we can find all the e's if we know the E's, and from 
§ 19 (i.) that we can find all the A's if we know the A's. It follows from (40) that in 
this latter case we can at once obtain the e's in terms of the o-'s. 

(iii.) We use the notation of § 7, and we also write 

— OfjE: coefficient of Sj (as auxiliary) in (e^_ 1? 
so that 

/ti = Oif/>i, . (74) 

e jtj .= i. (75) 

It should be observed that 9 ftj is not equal to 6 jtf . The 0's may be known directly, 
or, as is shown in (83), we may be able to obtain them from certain of the A's. 

18. Formula of Progression. — The quantities which we want to find are the 

improved values 

of cf , using S l9 § 2 , S S} . . . , 

of S x and S , using $ 2 > S Si S^ ... , 

of 4? ^u an( i § , using S 3i c! 4 , S 6} ■. . . , 
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and so on. There is a formula connecting these, which makes it unnecessary to deal 

with more than the first quantity in each row; or, if we deal with them all, the 

formula can be used for checking the results. (An example is given at the end of 

§ 15 of " Reduction/') 

We nave 

(e f )j-i = 8f-Of t j8j+ terms in $ j+1 , § j+2} ... . 

By (XV.), this is the improved value of S f —9 f jS p using S'$ after <^ ; and therefore, 

Dy l^iLllI.j, 

( e /)i-i = ( € f)j-~ e f,i( € j)i = ( e /)j- /,j E i- • ...... (76) 

Re-arranging, and replacing/ by j—l 9 /— 2, .../+1, and remembering that, by (75), 
ff = i 9 we have 

( e /)i ~( e /)j-i ~ fy<j % 



( € /)/+l""( e /)/ ™ ^/,/+lE/+l, 



Hence, by addition, 



( € f)j - a/^/ +0 a/+^ . . . . . . (77) 

19. Mean Products of Error (Alternative Formula), — (i.) By (77) and (38), 

(^f, g )j '= m.p.e. of <^ and (e^) 7 - 

= m.p.e. of Og and 0/yJiy+6U/ +1 Joy + i.+ ... +0/jJtL* 

~ /,/(^/)/+ e /,/+i(V/+i)/+i + ••• +®Aj(\>j)j 

The summation has to be made from £ =fto t = /. But, if <7>/, we Bee from (37) 
that it is sufficient to make the summation from t ~ g. Hence, using c ' t — f g " to 
denote summation from £ =/ or from t = g, according &s f or g is the greater, 
we have 

(fy.^/j = ^ @f,t(\t)t- • • • (78) 

t=ff9 , 

But, by putting # =y in (78) (or taking the m.p.e. of ^ and each member of (76)) 
and then replacing f and / by g and £, 

\ "-g*t)t Vg^l\.p . . . . . . . . . . y§ \y J 

Hence, substituting in (78), 

V^f'ph = ®f,t@g r tAp * • T * ♦ t (80) 
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If we can obtain the 0's and the A's in a simple form, we thus have a workable 
formula for calculating the A ? s, and thence, by (40), for determining the es in terms 
of the <x s. 

(ii.) From (80), using (II.) of §6, we get the m.ss.e. and m.pp.e. of the improved 
values of any l.cc. of the 8 J s. Let 

w = b § + Mi + . . . + b t S l9 w f = c ^ + c^, + . . . + c^, 

and let the improved values of w and w\ using <fs after S j9 be # and as'. Then 

m.p.e. of x and x f = 2 (Mo,« + Mu+ •••) (co#(u + c i$m+ •••) ^ 

- 2 (Mo.«'+Mi,«+ ...+Mm) Mm+^u+ — +c«e M ) A » ( 81 ) 

m.s.e. o£ x = 2 (Mo,{ + Mu+ ••* +&A t) 2 ^- . . . . . (82) 

i = o ' ■ 

(hi.) We have assumed that the ? s are known. If they are not known directly ? 
but the values of (\f tt )t are known, then, by (79), 

Substituting in (80), 

(^fig)i == ^ v^fit/t \\,t/tl^t (84) 



Also (77) is replaced by 






Application to Self-Conjugate Set. 

20. Preliminary. — (i.) We have now to apply the preceding results to the case in 
which the us are a self-conjugate set, so that (u r ; u s ) = (r^s), (u r ; u r ) = I, y r = u r . 
We take the § 7 s to be successive differences of the us, commencing with a difference 
of order 0. The <5's to be used as auxiliaries will be those following Sj ; the ( ) 7 
will usually be omitted. We shall take the number of us or of <5s to be m, so that 

m = Z+ 1. 

(ii.) There will be three cases to be considered; advancing differences, and central 
differences for m odd and for m even. For advancing differences the us will be 
taken to be u , u u ...u m __ 1 . For central differences we shall write m = 2^ + 1 or 
m = 2n ; and the ?/s will be u_ n , u_ nfU ...u n and u__ n+1 , u_ n ^ -.. w w respectively. We 
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shall require the following m.pp.e., which can be obtained from ordinary difference 
formulae. 

m.p.e. of A f u and A% = ( — ) f ^ g {f+g>f) 9 . (86) 



5? 



5J 



?> 



J) 



?> 



J) 



^u a „ <f% =(-) f "(2f+2g ) f+g) ) (87) 

o V/q ,, fxo %Iq =: U, (ooj 

M^"X „ M^-X = (-) / ^(2/+2^~2 ? /+^~l)/(2/+2^), . (89) 

3"-Hj» „ ^"S =(-) / ^(^+?flr-2,/+ fl r-l), .... (90) 

P-% „ m^" 2 ^=0, ' . (91) 

M^S „ M^" 2 ^ = (~) / ^(2/+2^^4 J /+^~2)/(2/+2^-2). (92) 



(iii.) For advancing differences we shall have 

S f = A% , ey = AA; . 

The formulae will be marked (A). 

(iv.) For central differences the two cases of m odd and m even must be considered 
separately ; but it will be found that, when the formulae relating to v 09 §\ 9 . . . 
(m odd) and to §v h S% 9 ... (m even) are properly expressed, they are practically 
identical in form, as also are those relating to julSv 09 ju$\ ... (m odd) and to 
fA.v i9 fitiPvi, ... (m even) ; and the latter correspond to the former with certain inter- 
changes of ( ] and [ ). We therefore, for ^§v ih ^\, ... and w h ^ 2 v h ... , replace 
6, E, A, X, by <f>, I, M, jul 9 with the appropriate suffixes. 

(v.) For m = 2n+l it will be seen from (88), taken with (XIV.) of §6, that the 
differences of even and of odd order can be treated independently. The $& will be 
u 09 S 2 u 09 . . . § 2n u in the one case and julSu , /ul§*u 09 . . . julS^^Uq in the other. We shall 
denote these by § 09 $ 29 ... § 2n and S l9 <J 3 , ... ^ ml respectively, and shall take j to be 
2k or 2k + 1 for the former and 2&-1 or 2k for the latter. The subscripts of the 0's 
and the ^s will' be modified accordingly ; i.e., 6 2fi2k will mean the coefficient of -S 2k in 
( € 2f)2k-2> a nd similarly for 2 /-i,2£~i- The formulae for the two cases will be marked 
(B) and (C) respectively. 

(vi.) Similarly for m = 2n we see from (91) that differences of odd and of even 
order can be treated separately. The $s are S l9 <5 3 , ...S 2n ^ in the one case, and 
S 09 S 29 ... S 2n _ 2 in the other, where S 2f _ x = SV- 1 ^ S 2f __ 2 = ^~% ; and/ is taken to be 
2&-1 or 2k for the former and 2k-2 or 2k-l for the latter. Also ft^,^ means 
the coefficient of -^.„ 1 in (<? 2/ _ 3 )aft-3 ; and similarly for <£ 2/ „ 2 , 2 *_ 2 . The formulae will be 
marked (D) and (E). 

vol. ccxxi. — A, 2 I 



224 DB. W. F. SHEPPAED ON 

(vii.) Writing 

Win R • n yfr \ = 1 4- * * ^ •-- *- »(«+l)'^(^+l)"> . 

p.\fr... p [p-h 1 J . ^ (^ + I j ... 
where a is a negative integer, it should be noted that 

2P{_n,/3;l,*} = fc&^, (93) 

and that, if 

then # 

l>, n] Lx> n] LV» W J Lx. raj . 

(viii.) For the central-difference formulae it will be convenient to write, if r and s 
are both even or both odd and s ;> r, 

{ , >r }=(X).(«±^ «)(».*» + **•), ,. . . . . . . (95) 

T -j- 1 



{r, 4 = (l) r (r+l)( i r ' + i S;r) ;- • • (96) 



so that, if k^f, 



s 2 k+i, 2/+ 1 } - [/+& *](&>/) = [/+-M+i](*+i»/+i) > . , ( 97 ) 

s 9 j, o/i _ C/+i *](*»/) _ l f+h k-l](k-l,f-l) (qR) 

and, if Jc <; s, 

lJ * +1 '^ 41/ ~(2A+2)(«, A) (2s+l) (*,*)' • • • • ^^ 

1 (2&+l)(s, &) 

(ix.) The successive steps are as follows. The formulae for the es (the improved 
values of the differences) in terms of the differences have already been found in 
" Keduction" and "Fitting" ; they depend on certain theorems as to the coefficients 
when l.cc. of moments or sums are expressed in terms of differences. From these 
formulae we get the 0's, and also the E's ; and thence we get, in each case, the 
progression formula supplied by (77). This formula is not really necessary, but it is 
useful for checking. The A J s, i.e. the m.ss.e. of the E's, are found from (86)-(92), 
using (IV.) of §6 ; this results in certain hypergeometric series, to which we apply 
(93) and (94). We then get expressions for the A's, by (80). From these, by (40), we 

* 'Proceedings of the London Mathematical Society/ 2nd series, x., 474, 
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have the values of the es in terms of the ars ; and also, by § 10 (iii.), the values in 
terms of the us (which are identical with the ys). This completes the investigation 
of the improved values ; but we also want to see the extent of the improvement. A 
further section therefore gives the ratio of the m.s.e. of the improved value to the 
m.s.e. of the original value, in a form convenient for calculation. 

(x.) The m.p.e. of A f u r and A f u r+t is ( — )* (2/, f+t). Hence, in finding the general 
solution of our problem for the case of a self-conjugate set, we are also finding it for 
the case of a set in which the m.p.e. of u r and u 8 is of the form ( — ) sr (2f 9 f-t-s--r) 9 f 
being some positive integer ; for we can treat these us as the f th differences of 
members of a self-conjugate set. 

21. Improved Values. — (i.) Adapting §11 (hi.) to the case in which the ys are 
identical with the us, we see that, if w is any I.e. of the u'&, its improved value, 
using differences of order exceeding j\ is of the form Hp r u r , where p r is a polynomial 
of degree j in r. The improved values of the us and their differences are obtained 
from this by means of certain formulae, given in §§6 and 7 (iv.) of " Factorial 
Moments," for the expression of ^p r u r in terms of differences. The results are given 
in (12), (18), (19), (26), and (25) of " Fitting." From the first of these, replacing 
m+lbym, we have (f = 0, 1, 2, '...j) 

( A) ^-,;>/)o--'.y-/) ^{:^^;:;) A^. i ..... do.) 

and from the other four, altering k to k— 1 in the last, we have (f = 0, 1, 2, ... k in 
(102) and 1, 2, ... k in (103)-(105)) 

(B) Wv -'27*+, k + f)(k , 7. f) (2k + 2f+l, 2k) [n+j, 2s + l) #, . . 

{B} lVo-2 o V + s,k+J){k-3,k J)l 2k + 2s + 1> 2k)ln+h2f+l) SU( >- (1 ° 2) 

-^n^^ i^s^ x. . (1 o 2A) 

. . . (103) 
^ Ua+{ ' .SUi*-f2k+l\{2k,28}iim,2f]'* U °' ■ ' (103A) 

(d) sv-\= s i n (k+s-i, k+/-i)(k-s, k-f) j^ + y~?' **-}\ b i s i **-%■ 



S= 1 



(2k+2s-l, 2k-l)[n, 2/) 

. . . (104) 



fff-hu. + ( _ Y-f S i" *=* 2/+ 1 {2^, 2/} [jm, 2s) », . , v 

i+[ } sXis-f2k+l{2k,2s}[Jim,2f) d U *> ' .' ' (1 ° 4A) 



2 i 2 
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(E) ^ Vi - T<* +s - 2 , i+/ - 2 ) (*- s , *-/) ^±i^=|) (^Ij ^ 



% 



8 = 1 



. . . (105) 



. s = k+is-f2k{2k-l,2s-l}(f>n,2f—l] 5 



(ii.) From (i.) we obtain 



(A) 



(B) 



(C) 



(D) 



e 



f.j 



u 



2f, 2k 



02/- 1, 



2k -1 



e 



2/-~l,2£-l 



(E) 



$2/- 2. 



24-2 



. . . (105a) 

' a/; ( 2? ;y-i) K/+i)- • ■ • ■ l 106 > 



vt., {2&, 2/} [jm, 2fe+l) 



) 



{2k, 2Jc} [±m, 2f+l) 



Y-f (2^-1, 2/-1} (fan, 2ft] 
' {2h-l, 2ft- 1} (|m, 2/]' 



_ }*-/ {2A:-1, 2/-1} [j-m, 2ft) 
j {2ft-l, 2ft-l} [|m, 2/)' 

_\*-/ (2^-2, 2/- 2} (|m, 2^-l"1 
7 ~{2ft-2, 2ft-2} (im, 2/-1" * 



(iii.) Also, by putting/ = y in (101) and/ = ft in (102)-(105), 



(A) 



(B) 



(C) 



(D) 



(E) 



E,= A^ =' m 2\8j)W +1 'J\ ( m '* +1 } * 



% 



o> 



XT 

-Hit 



-^2* 



#*» 



2 



:f {2ft + l,2&+l} fa 28 + 1) 



= k {2&+1, 28+1} [i|-m, 2ife+l) 



<f%, 



0? 






2£-l 



U 



if i 2k > 2 ^> te 2 *1 #. 



= jb {2&, 28 } (|m, 2& 



^o? 






8 = *{2Ms} [£m, 2£) iJ 



2;i ._ 2 = M o J * % = S 4—}- ' T/i^ — y 4 MO Ui. 



= k {2k-l, 2s -1} (fan, 2ft- 1] 



(107) 



. (108) 



. . (109) 



(110) 



(111) 



(112) 



(113) 



(114) 



(115) 



(iv.) It has been pointed out in § 11 (iv.) that the differences of order j all have 
the same improved value. It follows that (112)— (115) are particular cases of 
(ill), expressed in terms of central differences ; the proper values being taken for 
m and for/, and u a being altered to u_ n for (112) and (113) and to u_ n+1 for (114) and 
(115). We can verify this by expressing the E's in terms of the differences of 
order y. For (A) we have 



AX = {(1 + A) — I}* - -' A% = A% s _;-(s~y, 1) A% s ,;, + ... . 
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Substituting in (ill), and rearranging the terms, it will be found that the coefficient 
of A j u r is 

{j+T ^ j) ^^M^M Fi-fn+j+r+lj + r + l ; l,2j + r + 2} 

= U + r,j) {m-r-l, j)/(m,2j +l~\; 
so that 

< A > E < = "!> + ^> j™m4 Ah> < U6 > 

Similarly from ( 1 1 2 )-( 1 1 5 ) 

(B) E. = "IT* ( n + k+r ; 2k ]i n +^ 2 *) ^u r , (117) 

r --n+]c (^ ? 4#+lJ 

(C) i._ ]= r T'fc±*±!^=^ (118) 



r = 



(D) E. x = 1" fe±*±T-l. 2^-l)(n+^r z ^j^ I ll^_ (U9) 

r«-»+* (m, 4A:— 1J 

(E) ]^ = Tig±*i <-l. 2t-2)(n + *-r-a > 2*-2) (j*-^.^-^). . ( 120 ) 

The identity of (117)-(120) with (116), the u's being altered as explained above, is 

easily verified. 

(v.) The formula of progression (77) takes simple forms if we attach the factors 

involving m to the <5's ; for m then disappears from the formula. 

(A) Writing 

A t = (m, « + 1) E f = (m, « + 1) ( A% )„ 

so that, in effect, we take S f to be (m y f+ 1) AA^, (77) gives 

(A) (m,/+l)A^ = % J (-y-f(t,/) ^f±=^ A t (121) 

t=f {Zt,t~L) 

For/ = 3, for instance, we should have 

(m, l)v = A —A r +%A 2 —%A 
(ot,2)Av = A^Az+iAs 

(m y 3)A\= A 2 -2A 9 

(m, 4) A 3 i; = .4 



3 



^ 



3 y 



where ^4 is the value of (m, 1) v for j = 0, ^4 X is the value of (m, 2) Av for / = 1, 
and so on. This may be verified by § 5 (i.)-(iii.) of " Pitting." 
(B) (C) Writing 

Pn = H«*, 2* + 1 ) E 2< , # 2< _! = (|wi, 2«] !«_!, 



20Q 



DE. W. F. SHEPPAED ON 



we have 
(B) 



(C) 



[tm,2f+l)& 



v, 



T ( \t-f {%h 2«/ J P /'199^ 

K J {2t,2t} u * ^ ' 



(¥">> 2/] p& f -\ - 



* = / 



«-/ j {2*-l, 2«-l} V *- 1 " ' 



(123) 



The first of these has been given, for /= 0, in "Reduction," § 15 (v.), p. 362 ; the 
notation of P it differing, however, in a factor Jm. 
(D) (E) Writing 

we ha^e 
(D) 



2*-l - 



[fan, 2t) E 2 ,„ 1? Q^ 3 = (im 9 2«-l] I 

t = 
t 



2t~2) 



/T7l\ 



{^2f)^h h =l\~r f ^^-~^P,-,, (124) 

(^,2/-!]^-^= 2* (-)'-/ ^=||^ ^_, (125) 



(vi.) The A J s are obtained from the E f s by means of (IV.) of § 6 ; e.g., for advancing 
differences, m.s.e. of Ej = m.p.e. of E> and A% . We can use either the values of 
the E's given by (ill)— (115) or those given by (116)-(120); for the former we 
require the m.pp.e. given in (86), (87), (89), (90), (92), and for the latter the value 
of the m.p.e. of two differences of order f as given in §20 (x.). Using the former 
method, we get the following results :- — • 



(A) 



A, 



s =». m- 



1 ( _ y-J ( s j) M+UL (m,s + l) , ■ -x 

1 ' lJJ ( S +i+lj-)(m, j + l) lA+Wj 



(B) 



A 



2k 



(O) NLm-i 



(2y\j)F{—m+j+l 9 2j+l; 1, 2/+2} 

(2j 9 j) (m-/-l) ! (2/+1) \/{m+j) I 

(2j 9 j)f(m, 2/+1]; ........... 

. = * {2/c+l, 2s+l} [|m, 2&+1) v 7 
(4A, 2*) J P T {-w + A,w + A+l, 24+i; 1, 2ifc+l, 2i+f} 
(4&, 2k)/(m, 4& + 1] ; 

'if (_y-* {2^, 2&} (j-m, 2s] (2s + 2k— 2, s + &— l ) 
• * ' {2A, 2s} (£m, 2k] 2s + 2k 



(126) 



(127) 



s = Tc 



\D) &2k-l 



{4.k-2, 2k-l)/(&fy . F {-~n + Jc > n + k+ l,2A-£; 1,2*+ 1,2*+ J} 
(4*— 2, 2*— l)/(m, 4*— 1] ; 

S" (") S - A HHrJ fl 2 ^ (2* + 2*-2, * + *-l) 

s = s 7 {2&, 2s} [%m,2k) x ' 

(4k-2,2k-l)F{-n + k,n + k,2k-i; l,2k,2k+%} 

^4/C ,&, -qaJ — jLjjiTyb) 4a7"~" 1 J J • ' * . . » . « . . . » . . 



(128) 



(129) 



REDUCTION OF ERROR BY LINEAR COMPOUNDING. 



229 



(E) M 



2k -2 



, = * ' {2&-1, 2s-l} (|m, 2&-1] 2s + 2£-2 
== (4*-4, 2£-2)/(4*-2) . #{-n + A, n + *, 2A-f ; 1, 2*, 2*-£} 
= (4&— 4, 2Jc~-2)f(m> £k — 3] 



(130) 



It would, of course, in view of the identity of the E's as stated in (iv.) above, have 
been sufficient to obtain (126) and deduce (127)-(130) from it. 

(vii.) The A's having been found, the \'s, i.e., the m.pp.e. of the improved values of 
A f u and A% , etc., are obtained by (80). The 6 7 s and $>'s being as in (ii.) above, 



(A) 



<=i 



X /^ = S e f,t%,tfa, t)/{m 9 2* + l], 

t=f,g 



» • 



» « « « 



(131) 



(B) 



t=s 



X s /.s,= 2 0a/,»0*.*(4*, 2*)/(m, 44 + 1], 



. (132) 



(C) 



P) 



(E) 



i=fc 



2 A. 



V-2f~l,2g-l ~~ *-> y>2f-~l,2t~l T'<ig 

t=f,g 



**-i.*-i(4*-2, 2*-l)/K 4*-l], . . (133) 



£ = £ 



V-i,2,-i = 2 2/-i,2«-i ^-1,2^1 (4*-2, 2*-l)/(wi, 4«-l], . . 



. (134) 



* = & 



^2f~2,2g^2 ~~ 2 < p2f-~2,2t~2 ( j>2g~2,2t~2 

t=f,g 



(4*-4, 2*-2)/(m, 4* -3]. . . (135) 



To find the m.s.e. of the improved value of any I.e. of the differences, or the m.p.e. of 
two such Lea, we apply (II.) of § 6, as in § 19 (ii.). Thus, for m = 2n+ 1, 

m.p.e. of b v + 6 iy a^ + h^v + ... + & 3 *<^\ and c 'V + c 1 /x^ + c 3 (?V } + ... + c 2 ^ 2 \> 



) 2fV2f,2t\i ^ C 2g&2g,2t 



2 <i 2 6 9/ 9/ «N"S CoA-o,K4«, 2*)/(m, 4*+l" 



f = o I / = 



+ 2 J , 2 
t= i 



2 6 2 /-i02/-i,2«-in 2 Vi^-i.2*-i ( 4i "" 2 > 2£-l)/(m, 44-1. 



(viii.) The A ? s having been found, the Vs are then given in terms of the cr's by (40). 
Using the expressions for the o-'s given in (26) and (3l)-(35), we get the following :— 



(A) 



A* 



=3 



V, 



2j Kf g 2j Uq . (lob) 



= 



L\ 2 
r/ ■--= 



(-y\ f ,^h V , sm« 



t—m 



£=0 



• • 



(B) 



(C) 



cf%, = 



Z<! 2 X 2/ , 2y /ior %+1 « ; arWj 



2 

= ° 



in* 



— fm 



. (136a) 



• (137) 



• . 



/X 



J 8 '- 1 ^ - 



= & 



£ = 3 to 



£ = — |wj 



f • » 



(138) 
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(D) 



(E) 



?f~ l Vi 



/x 



8*f-\ 



#=1 



- U = 1 J . 



« = 



£ = ^7» 



£ = — ^W 



§W» 



i = 



•|m 



. . . (139) 



(140) 



The expression in (136) differs slightly from that given for A f v in " Fitting," (15) and 
§ 5 ; see Appendix III. 

(ix.) Finally, we want to express the e's in terms of the us. This is done by means of 
the general theorem in § 10 (hi.), that the coefficients of the ys in any e are related 
to the m.pp.e. of this e and the e's in the same way that the v's are related to the e's. 
Thus we find that — 



(A) 



(B) 



r = m 



If A f v = 2 p r u ri then p 



r = 



r = 2 (r, g) \ f ; 

g - 



r = — n 



<7 = 



. (141) 



r = ??, g = k 

If ^ f v = 2 £> r ^ r , then p r = 2 [r, 2g)\ 2f>2g ; (142) 



(C) 



r ~ n 



g=Jc 



If haS^^Vq = 2 £V^ r , then p r = 2 (r, 2gr— 1]m 2 /.~i,2£~i ; . . . (143) 

= 1 



r = — »i 



P) 



r = »i 



</=& 



If <f/ fy = 2 jp r w n then p r = 2 [r-|-, 2gr-l)x a/ _ lf ^_ 1 ; . (144) 



?• = —71-+-1 



0=1 



(E) 



r = vi 



g~k 



If ^ 2/ \>* = 2 p r ^ r , then p r = 2 (r— J, 2gr — 2] fi 2 f-2,2g-* • • (145) 

= i 



r = — ?&+l 



For a comparison of these formulae with those given in " Fitting," see 
Appendix IV. 

22. Extent of Improvement (Central Differences).- — A question of practical 
importance is the extent to which the use of these formulae actually reduces the 
m.s.e. of some selected quantity, such as, for the cases marked (B), u or § 2 ^u . The 
m.ss.e. of the various improved values are found from (I3l)-(l35), by putting g =f. 
Comparing these with the m.ss.e. of the original values, for the central-difference 
formulae (which are the important ones for practical use), we obtain the following : — 



(B) 



m.s.e. of S 2 \ = 1 ' = k \{2t,2f) [fan, 2t + l) \ 2 (At, 2t) 
m.s.e. of <f% (4/, 2f) t =/ \{2t, 2t} [fan, 2/+ l)J (w, U + 1] 



1 ' = k U + l \ {2t, 2f} \ 2 



(m, 4/+l]*=/4/+l \{2f,2f}\ 

{ot 2 -(2/+1) 2 } {m 2 -(2/+3 ) 2 } ... (m 2 -(2*-l) 2 } . 
{m 2 -(2/+2) 2 } {m 2 - (2/-f~4) 2 } ... {m 2 -{2t) 2 } ~' 



(146) 
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(C\ m.s.e. of pff-ivo 4f % k [ {2£~-l, 2/-1} (fm, 2£]\ 2 (4*-2, 2t-l) 

1 j m.s.e. ofpf'- l Uo (4/-2, 2/-1) ,=, l{2«— 1, 2*-l} (|-m, 2/]J (m, 4*-l] 

J/ ' = * 4S-1 f {2t-l, 2/-l} \ 2 

"(m,4/-l].t,4/-lW-l,^-l}J 

{m 2,( 2/+1 )2 } {m »_( 2/ +3 )» } _ { m «-(2*-l)»} . (u7 v 

{m 2 - (2/) 2 } {m 2 - (2/+2) 2 } ... {m 3 - (2£-2) 3 } ' V ' 

(T)\ m.s.e. of S^\ 1 '£,* f {2t — 1, 2/-1} fa 2*) \ 2 (4*-2, 2*-l) 

1 j m.s.e. of^~S (V-2, 2/-1) *=/ 1{2*-1, 2£-l} [^m, 2/)J (m, 4*-l] 

1 ^ U-l \ {2t-l 9 -%f-l} \* 

(m 5 4/-l]^/4/-l|{2/-l,2/-l}j 

{m 2 -(2/) 2 } WM2/+ 2) 2 } ... |m 2 -(2*-2) 2 } . (u8) 
{m 2 -(2/+l) 2 } {m 2 -(2/+3) 2 } ... {m 2 -(2£-l) 2 r V ' 

tm m.s.e, of y.&f-\ 4/- 2 l = k f{2*-2, 2/-2} fe 2^-l] \ 2 (4*-4, 2*-2) 

j m.s.e. of ^ 2/ ~S (4/-4/2/-2)tt/l{2*-2, 2*-2}(^2/-l]j (m, 4*-3] 

= 4 /~ 2 v* 4 ^" 3 f {2^ — 2 ? 2/-2H 2 
(m,4/-3]^/4/-3l{^-2,2/-2}J 

{m 2 -(2/) 2 } {m 2 -(2/-f 2) 2 } ... {m 2 -(2£-2) 2 } , , * 
{ m 3_(2/-l) 2 } { m ' 2 -(2/+l) 2 } ... {m 2 -(2*~3) 3 } V ' 

23. Smoothing. — When we have^ a table containing a very large number of us 9 a 
common method of procedure is to use a formula involving a limited number of terms 
and to apply it to successive sets of the us for the purpose of obtaining a table to be 
substituted for the original table. Thus we might use a formula involving 2n + 1 
terms, and apply it to u 09 u l9 u 29 ... u 2n for finding a new value w n9 then to u u u 29 u 39 ... 
tt 2n+ i for finding a new value w n+l9 and so on. These values having been obtained, a 
differenced table would be formed ; but, as by hypothesis the true differences of order 
exceeding j are negligible, the table would only go up to differences of order j. There 
are two cases to be considered. 

(i.) If our object is to obtain as accurate values as possible for the w?& 9 consistently 
with our using only the specified number of u's for each, the most accurate values 
would be the u's given by the formulae considered in this and the preceding papers. 
It should, however, be observed that the differences of the w's are not the same as the 
A f v 0i § 2f v 09 etc., occurring in those formulae. Suppose, for instance, that we replace u 
by its improved value v obtained by means of the (B) formula involving u_ n9 u_ n + l9 
. . . u n9 and replace u x by the improved value v l obtained in a similar way. The 
resulting value of v x — v will involve the 2n + 2 lis from u_ n to u n+l ; but it will not be 
the same thing as the improved value v$ obtained by the (D) formula involving these 
us 9 and its m.s.e. will therefore be greater than that of the latter. 

(ii.) If our object is to obtain a smooth table of the w's as a whole, we could do this 
by obtaining as accurate values as possible for the differences of the w's of order/. 
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The formula which would have to be applied to the us in order to obtain this result 
can be constructed without difficulty. The important thing to notice is that, if we 
alter the differences of the us and then obtain the ws from the altered differences by 
summation, the resulting values must be such as can be legitimately substituted for 
the us. Suppose, for instance, that/ = 2&+1, and that we use 2n+l us for each w. 
The formula for w will have to be of the form 



w = u + c 2k+2 S 2k+2 u + c 2k+ ,S 2k+i u + . . . + c 2 J 2n u ; 



and this will give 



8*+hi>i = 3 8 * + S + c» + ^ + S+^ + 4^ + S+ - +c 2 J 2 » +2k +%. 

The problem of determining the cs so that the m.s.e. of ^ l w^ shall be a minimum 
is the same as that of determining the coefficients in the improved value of S 2k+1 u% for 
j -= 4&+1 or 4&-f 2, m being 2^-f2& + 2 ; and the solution of this problem is given in 
§ 21. Thus, in terms of sums, (139) gives 



§ 2k+1 w. 



* 



= 1 



*=s ~(»+A+l) 



The X's having been found, we shall then have, by summation. 



w, 







g = 2& + l 



Xy -s — 2j A 2 ^ +lt 2^_i^a" %q J (T 



ft+2, 



'Ui 



= 1 



$ = - (n + ft+1) 



The ratio of the m.s.e. of <f* +1 w$ to that of S 2k+l u h is given by (148). 



Appendix I. — The Correlation-Determinant. 
1. The m.p.e. of A and B being denoted by {A ; B), let 



e 



{A ;A) (A; B) {A ; C) . . . 
(B;A) (B;B) {B ; G) . . . 

(C;i) ((7; 5) (0;0) . . . 



We call this the correlation-determinant of A, B, C 9 — 

2. The elements of this determinant may be regarded as obtained as follows. We 
first take a representative collection of N A values of the error of A ; N A being usually 
indefinitely great. Then, for each of these values, we take a representative collection 
of N B values of the error of B ; the resulting N A collections will all be alike if the 
errors of A and of B are independent, but not if they are correlated. This gives 
N A N B combinations of an error of A and an error of B. For each of these we take 
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a representative collection of Rvalues of the error. of C; and so on. Thus finally 
we shall have N= N A N B N C . . . combinations of an error of A, an error of B; &c. 
Numbering these 1, 2, ... 'N, and denoting the errors of A, of P, of C, ... by 
a, &, c, ... , the combinations will be a u b l9 c l9 ... , a 2 , 6 2 , c 2 , ... , ...<%, 6 V , c^, ... ; and 
we shall have 

(A ; A) = (a 1 a +a 3 a + a 8 a +...+ <% 2 )/iV, 
{A ; B) = (a 1 6 1 + a 2 & 2 +a 3 6 3 +...+a iV 6 iV r)/iV r 5 

&c. 

3. Substituting these values in 0, we find that, if there are m of the quantities 

A. JlS. O, .... 



J -^J 



iV m e 



C^l + €b 2 + C£ 3 + . . . 



a 1 b 1 + a 2 b 2 -\-a s b 3 -{- ... a 1 c 1 + a 2 c 3 +a 3 c 3 +... &c. 
a 1 6 1 + a 2 6 2 + a s & 3 + ... 6^-f & 2 2 + & 3 3 +... &iO x + 6 2 c 2 4- & 3 c 3 + . . . &c. 

a x c x + a 2 c 2 + a 3 c 3 -f . . . b x c Y + b 2 c 2 + 6 3 c 3 + . . . c x 2 + c 2 a + c 3 2 + . . . &c. 



= (aj>>fi 3 . . -) 3 +. (^l^W • •)* + («l&3^4- • -) 2 + (^A^V . -) 2 + • • • 



where (o^Cg...) denotes 



hboK... 



1 V 2 V 3 



^1 ^2 ^3 • • • 



• • 



Hence is not = unless each of the 



determinants (a&c.) is = 0. This would be the case, for instance, if A were a 
constant, so that every a would be 0, or if there were a linear relation connecting 
the errors of A, P, 6 Y , .... 

4. Let <3> be the correlation-determinant of A, P, ... P, Q, ... , and "^ that of 
,4, P, ... P. 

(a) Suppose that ^r = 0. Then, by § 3 of this Appendix, each of the determinants 
(a& ... p), where a, 6, ... p are the errors of A, P, ... P, is 0. But these are the 
minors of the qs in the determinants (ab . . . pq) ; and therefore these latter 
determinants are 0. Proceeding in this way, we see that the determinants 
[ab...pq...) are all ; and therefore $ = 0. 

(6) Hence, if $ is not = 0, ^ is not — 0. 



Appendix II. — Frequency of Correlated Errors. 

1 . Let u 09 u ly u 2 , ...Ui and y , y l9 y 2y ... y x be two conjugate sets. Denote the errors 
of the us by O , Q u 29 ... Q l ; and let the resulting errors of the ys be <j> 0y <f> l9 <p 2y ... # z . 

2 K 2 
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Then, on the assumption of normal correlation of errors, the frequency of joint 
occurrence of these #'s is proportional to 



where P is a homogeneous quadratic functio~ ~~ .... 
'* x '" ^'s are the partial differential coeffici 



ction of the 0*s. We want to prove that 

(i.) the ? s are the partial differential coefficients of \P with regard to the & 
and conversely ; 



s, 



(ii.) P = e <t> + 6 1 <f> 1 + e 2 </> 2 +...+0 l <t> l ; 

(iii.) P = V r o.o^o a + 2^o f i^i + V r i.i i 3 +---+V r i,z^ a » where ^ is the m.p.e. of ^ and 
^ ; and similarly 

(iv.) P = 7T OjO ^ o 2 + 27r o>1 ^ o 0i + '7ri,i^i 2 +*-*+' 7:r z,^ 2 ' where ^^ is the m.p.e of u f and ^. 



2. Suppose that 



P = a OtO e o 2 + 2a Oil e o Oi + ^i t A 2 +...+a ltl l 2 i 



and let us, without making any assumption of conjugacy, write (/ = 0, 1, 2, ..J) 



Vf 



% 0^0 + a f s 1^1 + % 2^2 + • • • + a f % l U h 

error of y f 



2 



dP/de* 



Then, writing the subscripts in the order/, 0, 1, 2, ... I, 

P = a Af 6/+2a ftQ e f 6 + a Qt(i eo*+...+aiA 
= <pf/a f ,f+Q 9 

where $ does not contain fy. 

3. The mean value of <j> f 6 g is NjD, where 



N„ 



D 



<p f 9 g exp — JkP . c?6^ c£0 o c?^i • • * ^> 



J ... exp— ^P .d6 f d0 d6 1 ... c?0j, 



the integration being in each case from — oo to co . If we write 



then 



n„ = 



9 



*> J 



^ = 4>fl\/ a M> 

>/^ exp— J^ 2 . exp— J$ . d^ c£# cft^ . . . c£fy. 



(a) First, suppose that g is not =/. Then, integrating with regard to \fr t 



N g = 0. 
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(5) Next, suppose that g = /. Then 



N f = l Kf ■ 



<f)f y<py — CCf t q6 — Cb* i6i — ... — Cb* jOi) GXp — 2*JT . d Vf du dOi * . . d6 t 



= l/(%f,f. ... I #/ exp — JjP . d0yd0 d0i...d0i t 

by (a). Hence 

N f = ijs/a^f. ... \\f/ exp— ^\f/ . exp—^Q . d\p- d6 dB A . . . d6 t 



Also 



= \/ (2ir/a fif ) . ... exp — ^Q . c?0 o dfy . . . dO^ 



D = 1/a/%/ • 



exp-^.exp-W.^cWo^...^ 



Hence 



v 7 (27r/ay,y) . exp— f^ . c?0 o c?0x . . . <## z . 



tf/D = l. 



4. Hence the ys are related to the us in such a way that 

m.p.e. of y f and u g = (g?*f) or l(g =f); 

and therefore the w'g and the ?/s are conjugate sets ; which proves (i.). It follows 

that 

a ftg = m.p.e. of y f and y g . 

This proves (iii.); and (iv.) is the similar result which we should have obtained by- 
expressing P in terms of the ^'s. Also 



P 



^o,<A> + 2a o , 1 o 1 + a lll 1 +...+a Ll 6 l 
Qq \<h, o#o + a o, A + »o, 2^2 + • • • + ^0, fil) 



+ ... 



which proves (ii.). 






Appendix TIT. — Improved Advancing Differences in terms of Sums. 



The expression for A f v given by (136) differs from that given in (15) and §5 of 
" Fitting," in that it involves S"^, S" 3 ^, 2 //3 tt 2> ... , instead of S"^ S //2 w , S //3 tt € , — 
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The new expressions are more convenient for calculation and for tabulation, since the 
coefficients are rather smaller and are symmetrically placed about a diagonal. For 
j = 2, m = 13, for instance, the formulae given by " Fitting/' § 5 (ii.) 5 are 

100H' = +693S 1 -198S 2 + 22S- 4 , 
lOOlA^o = -23l5i + 885 a -ll^b 
1001A 2 i; = +35S 1 -15S 2 +2S 3 , 

where S 1 = 2X, &2 = 2"%>, # 3 = 2'X- If we write 2, = 2\ 2 2 = 2'%, 
2 3 = 2 //3 w 2 , these become (by (136), or by writing /Si = 2 l5 $ 2 = 2 X -f 2 2> 

A^g =:: 2^i "T* £2Li2 1 ^3/ 

1001v = +5172 1 -1542 2 + 222 8 , 
100iAv = -154S 1 + 662 a -llS 8 , 
1001AX = + 222 1 -112 2 + 22 8 . 
The symmetry of the coefficients is due to the fact that 

co. ?! tg + l u g in A f v = \ ftff = X^ = co. 2" /+1 i6 / in A% . 

For any particular value of m there will be only i{j+l) 0" + 2) coefficients to be 
tabulated, instead o£(j+l) 2 . 



Appendix IV. — Formula in terms of it's. 

(i.) Formulae for A f v , &c., in terms of the us have already been given in (15a), 
(21), (22), (29), and (28) of " Fitting " ; and the results in (I4l)~-(l45) of the present 
paper can be checked by comparing the different expressions for the coefficients of the 
us. We should require to use the following identities : — 

(r + h, h) - (r, 0) (h, 0) + (r, l) (h, l) + (r, 2) (h, 2)+ ... , 
(r, 2K] = (0, 2A] + [r, 2) (0, 2&-2] + [V, 4) (0, 2A-4J + ... , 
' [r, 2A-1) = (r, l] [±|, 2A-2) + (r, 3] [±J, 2h-4) + (r 9 5] [±£, 2A-6) + ... , 
(r-h 2fc-l] - [r-i 1) (0, 2A-2] + [r-i 3) (0, 2fe-4]-H>-£, 5) (0, 2A~6]+ ... , 
[r-h 2A-2) = [±i 2A-2) + (r-i 2] [±*, 2A-4) + (r-i 4] [±J, 2fc-6) + ... . 

(ii.) Taking, for instance, the formula for <5 2/ v when m = 2^+1, (21) of " Fitting " 
gives (replacing tbyf) 

(t>) -( )flm {2k+l,2f+l} %* ( Y {f+l)(h+l){2k+l ,2h+l}(r ,2fi] 



[fcw, 2/+1) ' fc = v ' /+A+i (£m, 2fc+l] 
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Hence we find that 



(P,k- (*>,k- 2 - 1 )2w»^+ 2 ) [ ^^ +1)h i^ J (^,2/^+1] ' 

\rrrn v /a ; = / \ ^/ j-^ 2/+1) ft = o (|m, 2A-+1] 

and therefore 

(^>o)» = (-)^im 'i* (2*+&) J 2t ' 2 {\ r Tl^^Vi 9 > 
\ i'wa \ / z 4 = f \ ^ [£m, 2/+l)(£m, 2# + l] 

where 

= i^{— £+0, S'+^+J, i ; 1, —fon+g-bl, %m+g+l} 

(jm, 2£+l][|m, 2gr + l)* 

This expression for ((^0)2* will be found to he equal to X 8/i a, as given by (132) of the 
present paper, so that the formula in " Fitting " agrees with (142). 



tmmmmmim 



